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We use a density-functional based tight-binding method to study diamond growth steps by 
depositing dicarbon species onto a hydrogen- free diamond (110) surface. Subsequent C2 molecules 
are deposited on an initially clean surface, in the vicinity of a growing adsorbate cluster, and finally, 
near vacancies just before completion of a full new monolayer. The preferred growth stages arise from 
C2n clusters in near ideal lattice positions forming zigzag chains running along the [110] direction 
parallel to the surface. The adsorption energies are consistently exothermic by 8-10 eV per C2, 
depending on the size of the cluster. The deposition barriers for these processes are in the range 
of 0.0-0.6 eV. For deposition sites above C2n clusters the adsorption energies are smaller by 3 eV, 
but diffusion to more stable positions is feasible. We also perform simulations of the diffusion of C2 
molecules on the surface in the vicinity of existing adsorbate clusters using an augmented Lagrangian 
penalty method. We find migration barriers in excess of 3 eV on the clean surface, and 0.6-1.0 eV 
on top of graphene-like adsorbates. The barrier heights and pathways indicate that the growth from 
gaseous dicarbons proceeds either by direct adsorption onto clean sites or after migration on top of 
the existing C2n chains. 
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I. INTRODUCTION 

The growth of ultra-uanocrystalline diamond films (UNCD) from C2 precursors produced by Cgo fragmentation in 
hydrogen-poor plasmastl has recently attracted attention because of high growth rates and resulting good mechanical 
and electronic properties of the films. „ PI-, 

There is evidence that the growth proceeds mainly on the (110) face.H Previous studiesOu explored growth by 
dicarbon on the hydrogen-terminated (110) face without hydrogen abstraction by way of insertion of C2 into C-H 
bonds on the surface. In the present work, we consider deposition steps onto the clean diamond (110) surface without 
hydrogen participation. Starting out from a clean surface we investigate the local atomic configuration arising from 
the adsorption of a C2 molecule. Subsequently, more C2 molecules are deposited in the vicinity of a previous adsorbate 
cluster. By comparing the total energy of these structures, we identify preferred growth stages as those arising from 
C2ri clusters in the form of zigzag chains running along the [110] direction parallel to the surface. 

We also perform simulations of the diffusion of Gp|-polecules on the surface in the vicinity of existing adsorbate 
clusters using a constraint minimization technique.l3 13 The barrier heights and pathways indicate that the growth 
from gaseous dicarbons proceeds either by direct adsorption onto clean sites or after migration above existing C2n 
chains. 

The computational method is briefly outlined in Sec. IL In Sec. Ill, the properties of clean diamond (110) surfaces 
are calculated using density-functional based tight-binding method and compared with more elaborate ab initio results. 
The initial adsorption steps for diamond growth are studied in Sec. IV, followed by an evaluation of diffusion barriers 
for C2 on diamond (110) in Sec. V. We then analyze some molecular dynamics trajectories in Sec. VI. Finally, we 
present conclusions in Sec. VII. 

II. COMPUTATIONAL METHOD 

In performing our atomistic simulations, we apply the density-functional based tight-binding method (DFTB).i In 
order to correctly take into account effects of surface polarization we have used the self-consistent charge extension 
of the method.El It has been applied successfully to diamond systems including studies of surfaces and diffusion 
problems.lijO It reproduces total energy differences between various carbon bulk structures within 50-100 meV/atom 
compared to ab initio calculations (Ref. HJlJ). For surface related effects, the errors are determined to be less than 
0.2 eV/atom. 

The surfaces are simulated using a 2-dimensional slab-geometry with varying thickness. The bottom layer is 
saturated by a fixed monolayer of pseudo- hydrogen atoms which do not mutually interact. The lateral extent of the 
cell varies from 3x3 unit cells up to 8 x 3 unit cells, where the bigger extent is along the chain direction. This 
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direction is assigned as a;-axis. For all atomic structure calculations, we used the F-point approximation to sample 
the Brillouin zone, which amounts to a fc-point sampling with as many points as real-space unit cells. We let the 
atoms relax in a conjugate gradient scheme. 

The diffusion and adsorption runs are done using a constrained conjugate gradient method. The^fosces during the 
conjugate gradient minimization are modified by method described by Ciccotti et al. and RyckacrtJ3n3 The center of 
mass of the C2 molecule is moved to a given direction with steps of 0.1 A. Because the constraint is to the center of 
mass movement, the C2 may rotate or dissociate freely. 

III. CLEAN DIAMOND (110) SURFACE 

The clean (110) surface of diamond as obtained from bulk cleavage has one dangling bond per atom pointing at 
an angle of 19.5° away from the surface normal. The atoms are arranged in zigzag chains flat on the surface directed 
along the [110] direction. 

For initial studies we used six bulk monolayers of carbon. We find that only negligible relaxations take place in the 
lower two layers and therefore we use only four carbon monolayers for the remaining calculations. The applicability 
of this size of the super-cell is tested by comparing the minimum energy geometries of a C4 adsorbate on the (110) 
surface with super-cells consisting of four and six monolayers, respectively. The comparison between fully relaxed 
structures calculated with both number of monolayers yields very small differences in the atomic positions of about 
0.015 A in both the lateral and vertical positions of all corresponding atoms, including those of the adsorbate,— 

There has been some controversy in the literature about whether the C(llO) surface reconstructs or not.E^O We 
find here in agreement with ab initio studiesEj and experimentstS a symmetric flat (1x1) surface. The upper surface 
layer moves inward by Azi — —0.16 A, and the next layer outward by Az2 = 0.02 A. The chains in the top layer are 
slightly straightened by a relaxation of 0.1 A towards the chain axis, resulting in a bond angle of 122.0° and a bond 
length of 1.44 A. This geometry is in excellent agreement with the ab initio results in Ref . |l^, summarized in Table 1. 
Taking into account the 0.015 A offset in the equilibrium bulk bond length do, all bond lengths agree within 0.01 A, 
or 0.6% of dn- 



IV. ADSORPTION AND ENERGETICS OF SMALL CARBON CLUSTERS ON (110) DIAMOND 

The energetically most favorable cluster configurations after repeated C2 additions are summarized in Fig. ^ up to 
Cs on (110). In the following, a more detailed description of the surface cluster geometries will be given. 

A. Initial C2 deposition 

In order to sample the energy landscape above the clean (110) diamond surface for C2 adsorption we place a C2 
molecule in vertical orientation near the surface on a hexagonal-like set of points above the atoms and bond centers 
of the two topmost monolayers. By symmetry, only seven positions remain unique. The lower atom of the molecule 
is placed about 2 A away from the nearest surface atom. A conjugate gradient relaxation from each of the lateral 
starting positions shows that the molecule is either refiected from or adsorbed onto the surface. The refiections occur 
for starting positions directly above the atoms and bonds of the top monolayer. From all of the starting positions 
not directly above a top-layer chain the C2 molecule is bonded and forms a bridge above the "trough" between two 
adjacent top-layer chains. The deposition proceeds in two stages, see Fig. ||. 

Initially, the C2 sticks with one end to the nearest surface atom with an inclination of about 45° to the surface 
normal. There is a very low energy barrier of order 0.1 eV towards the final adsorption stage in which the molecule 
bonds symmetrically in an orientation corresponding to the diamond-lattice. At this final stage, both adsorbate atoms 
are 1.0 A above the top monolayer and are 3-fold coordinated with a bond length of 1.38 A between them. Each 
adsorbate atom forms 2 bonds towards the surface: one bond corresponding to the diamond lattice with a length of 
1.49 A, the other being a stretched bond of 1.91 A towards the adjacent atom in the same surface chain, see Fig. |l| 
(a). The result is the formation of two 5-fold rings with a common bond formed by the adsorbed molecule. There is 
a slight lateral pinch contraction of the surface monolayer. All its atoms remain bonded to the subsurface, which is 
indented below the adsorbate by about 0.05 A. 

The binding energies of the C2 molecule to the clean (110) surface are fairly high, as can be seen in the fifth column 
of Tab. ||. The energy gain is mainly the result of forming four bonds to the surface which yields about 2 eV for 
each bond, which is plausible by its similarity to the atomic binding energy of 2.3 eV/bond as obtained by the DFTB 
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method for bulk diamond. Breaking the initial triple bond within the C2 molecule and the stretch of surface bonds 
near the adsorbate offsets the result to the adsorption energy of 8.1 eV, listed in the table. 



B. Addition of C2 

As a next step, we have studied the effect of adding another C2 molecule near an existing C2 adsorbate. Obviously, 
the deposition onto a site more than one surface lattice spacing away from the initial adsorbate results in two isolated 
clusters of similar configuration as discussed above, unless a topological mismatch prevents the completion of the ring 
formation for a near site, as shown in Fig. ^ (a). In this case, the strain field introduced by the second adsorbate 
results in a bond switch for one of the backbonds, with an accordingly high adsorption energy of —7.2 eV. For the 
deposition onto a site of the neighboring chain, the geometry and energy is essentially the same as for the isolated 
case, see Fig. ^ (b). We note that in this case, the low barriers found in the initial adsorption process are no longer 
present, likely due to the small local strain field induced. This effect supports a rapid spread of such C2 adsorption 
sites across the surface. 

C. C4 clusters 

The highest gain in energy for the second C2 deposition is obtained at a site directly above to the first C2 molecule 
(see Figs. H (b) and The second molecule bonds at the neighboring diamond-lattice site along the [110] valley, in 
a diamond-like configuration next to the first one and forms a 4 atom long zigzag chain which amounts to a seed for 
the next monolayer. At the ends of the new chain, 5-fold rings are formed similar to the ones in the C2 case. The 
ridge of the adsorbate is a z-sliaped symmetric chain of 3 bonds with lengths of 1.40 A at the ends and 1.47 A in the 
center. The central two atoms are raised 1.3 A above the top monolayer, which corresponds to a slightly outward 
relaxation with respect to the ideal lattice sites. The end atoms are 0.3 A closer to the surface. The center bond 
is the common side of two adjacent 6-fold rings connecting the new chain to the surface monolayer. Topologically, 
the rings supporting the C4 adsorbate above the surface form a pyracylene structure, which is the basic structural 
element of a Ceo fuUerene. A local pinch contraction of the surface monolayer occurs similar to the previous case. 
However, the outer bonds of the adsorbate are here with 1.6 A more closer in length to actual bonds, at the expense 
of the transition to the uncovered parts of the surface, which are now more stretched to a distance of 1.9-2.0 A. 
Furthermore, we observe the breaking of backbonds in the middle of the aggregate, resulting in a cluster of sp^-like 
coordinated atoms arranged in a dome-like configuration. The similarity to such a highly stable configuration as in 
a fuUerene explains the fact that this C4 cluster represents the highest energy gain for an approaching C2 molecule 
among the structures considered in Tab. ^ 

D. Cg and Cg clusters 

The next C2 adsorption to the C4 pyracylene-like adsorbate results in a Cg adsorbate, shown in Fig.|^ (c), with an 
adsorption energy of 9.6 eV. Further C2 adsorption yields an adsorption energy of 8.8 eV and a surface-Cs cluster, 
shown in Fig.^ (d). Thus, the energy gain is at least 8 eV for the repeated C2 surface chain addition, as shown up to 
Cg on the surface. We expect the adsorption gain to level off at about 8.5 eV for longer chains. 

The Ce adsorbate has four six-fold rings in the middle and two five-fold rings at both ends. The Cs adsorbate is 
basically identical to it, but has of course two more 6-fold rings along its length. As in the case of the C4 adsorbate, 
the underlying substrate atoms are raised above the surface and flattened to sp^-like coordination together with the 
adsorbed atoms. 

E. Adsorption barriers 

The energy barriers for the C2 adsorption (Tab. fourth column) are either zero or very low. The highest energy 
barrier in the case of C2 addition at the end of Cs is probably a finite size effect. 

The initial C2 addition at a diamond site in a "valley" on the surface makes the bonds shorter in the neighboring 
"valley" . This makes the adsorption of an additional C2 easier to the neighboring 'valley' at a diamond site, as can 
be seen by comparing the energies in third and fourth rows in Tab. |l|. 
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F. Surface defect formation 



There is a meta-stable energy minimum when adsorbing a C2 molecule on the top of C4 or Cg adsorbates (rows 6 and 
8 in Tab. y). In the meta-stable state {C2 + C2n) one of the C2 atoms is singly bonded to the surface C2Ti-complex, 
see Fig. | (b). The bonded C2 is only 2.0 eV lower in energy than free C2. We beheve this meta-stable minimum 
configuration plays a key role in the growth of (110) diamond. It enables the diffusion of the C2 molecule to the end of 
an existing C2n-complex growing diamond. If the C2 molecule is forced deeper onto a Cg, there is another meta-stable 
minimum energy structure consisting of one 7-ring, one 6-ring and six 5-rings (row 9 in Tab. ||). The energy gain 
from the gaseous C2 to this non-diamond-growth favoring configuration is 4.9 eV and the energy barrier towards the 
final meta-stable minimum is 1.8 eV. We believe that similar meta-stable defect structures form, when a C2 adsorbs 
to C2n with too high kinetic energy, approximately Ex > 3-5 eV, taking into account kinetic contributions to the 
energy barriers. 



G. Surface vacancy filling 

When the growth proceeds further, different C2n clusters along the same (110) surface trough will eventually meet 
and coalesce. Given that growth proceeds by C2 addition, the critical stage is reached just before coalescence, when 
there will be a gap between two cluster-ends corresponding to either 3 or 2 missing atoms. Assuming the clusters are 
seeded at random sites, both cases have equal probability, but quite different energetics for subsequent C2 additions. 

In either case we see the approaching C2 first in a meta-stable bridging configuration from a chain end to a bare 
surface site, and directly between the two chain ends, for the 3 and 2 site wide gap, respectively, see Fig. ^ (a) and 
(c). Both added atoms remain just 2-fold coordinated. After overcoming a barrier of 0.6 eV and 0.3 eV, respectively, 
the adsorbate extends the existing C2„ chain by another atom pair, see Fig. || (b) and (d). 

For the original 3-atom gap, a single-atom vacancy remains next to a still just 2-fold coordinated atom. Accordingly, 
the gain in adsorption energy is with 6.3 eV rather low. However, the remaining single-atom gap remains reactive, 
and may be filled at a later stage. 

The 2-atom gap yields a much higher adsorption energy of 10.2 eV, comparable to the high gains found in the initial 
adsorption stages. The final stable configuration is a contiguous chain with broken backbonds for atoms on either 
side of the top ridge. This structure is a bent graphene sheet with a bending radius of about 3 A. Before discussing 
its properties in the next chapter, we briefly sketch the other variants for surface vacancy filling. 

For the final stage of surface coalescence, we considered a nearly complete surface monolayer, with up to 3 consecutive 
atoms along the [110] direction removed. The filling of these surface vacancies results in energy gains between 6.9 
and 10.4 eV, as listed in the last rows of Tab. The filling of the last single-atom vacancy can take place either by 
by a single C atom adsorption to the vacancy, with an energy gain of 10.4 eV without a barrier, or by a C2 addition 
process, which has a rather high barrier from a meta-stable minimum at a gain of 6.8 eV to its completion at 8.4 eV. 
Furthermore, it leaves a singly bonded C atom on otherwise perfect (110) surface. The energy required to desorb the 
extra C atom is of order 8 eV. 



H. Graphitisation and rebonding 

As shown in the preceding paragraphs, the growth process by the C2 chain addition mechanism will eventually lead 
to coalescing chains, with broken backbonds on either side. In order to investigate the consequences of the broken 
backbonds for the surface stability during growth we have relaxed a model in which every other trough along [110] 
was covered with a contiguous chain, resulting in a C(110):(2 x 1) reconstruction, shown in Fig. ^ (a). The relaxed 
structure shows multiple bent graphene sheets along [110]. The bending orientation is that of a carbon nanotube of 
the (n, n) type, known as armchair tube.EZI 

It is knownJJaat graphitisation on clean diamond surfaces, namely on (111) and near (111) twin boundaries, leads to 
delamination.E3 We find for the present configuration that it is stable and does not debond. Furthermore, continued 
adsorption of C2 in the valley between two arches is possible without a barrier and more importantly, it causes 
the sp^-like atoms near the adsorbate to return to an sp^ configuration and rebond in the diamond structure, as is 
illustrated in Fig. (b). 
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V. SURFACE DIFFUSION OF Ca 



In order to estimate the influence of surface diffusion of C2 on the growth mechanism, we investigated some diffusion 



paths, as summarized in Tab. III. The associated diffusion barriers along the various C2-related diffusion paths are 



shown in the second column of this table, and the gain in energy in the last column. 



Generally, on the clean diamond (110) surface, the diffusion barriers are rather high, and exceed 3 eV (Tab. Ill 



rows 1-3). This is easily understood by the strong covalent bond that is formed between a C2 ad-species and the 
surface once the molecule reaches the surface. The only exception to such high barriers are for sites above existing 
adsorbates where the binding energy for further ad-species is low to begin with. C2 has the lowest diffusion barriers 
when it starts diffusion on top of an existing C2n complex. In this case, the barrier for diffusion along the adsorbate 
ridge is of the order of 1 eV. The C2 remains nearly vertically oriented, with one of its atoms bonded to one or two 
surface atoms throughout the diffusion path. The energy barriers are decreasing when the chain end is approached. 
The last energy barrier towards completing a chain addition step is only 0.6 eV. 

We note that the energy gain attainable for a C2 molecule by diffusion to the end of an existing C2n adsorbate 



is considerable, see Tab. [II, last column. These gains, when added to the adsorption energies found for the "top" 
deposition sites, as listed in Tab. ||, naturally result in the same total adsorption energies as those for the "end" sites. 
We have thus found two different growth channels, converging to the same growth mechanism. One is adsorption- 
dominated growth on nearly clean surfaces with deposition directly into diamond lattice sites, and the other diffusion- 
driven on surfaces covered densely with adsorption clusters. 

In the diffusion studies, we found that when the C2 is directed in both major diffusion directions, either along 
the chain direction of C2n, or perpendicular to it, it passes metastable minima before reaching the lowest energy 
configurations at the chain ends. Along the parallel diffusion path, the last metastable minimum consists of a C2 
fragment singly bonded to the surface, similar to the configuration shown in Fig. ^ (b). For the perpendicular direction 
away from the chain axis towards a neighboring surface valley the meta-stable minimum is such that the C2 atoms 
complete a 5-fold ring, with each one being singly bonded to the surface. 

While both intermediate and final diffusion barriers are incidentally the same for diffusion parallel and orthogonal 
to the C2n chain direction, the ultimate energy gain is higher by 0.6 eV for the growth-favoring parallel diffusion; this 
indicates high adsorption rates in either case, with a preference towards cystalline growth. However, the introduction 
of defects is quite easily possible under this regime, which helps to explain the rather small grain size found in the 
final material. 



VI. MOLECULAR DYNAMICS DEPOSITIONS 



Inspired by the adsorption and diffusion results discussed above we simulated the deposition of C2 on top of a Cg- 
complex on the surface using molecular dynamics directly, though within a rather short time span of just 0.1-0.5 ps. 
The goal of these runs was to investigate a diamond growth reaction as follows: 

Cgas . ^surf /-^surf / 1 \ 

2 "T ^2n ~^ ^2n+2- l-LJ 

We chose the initial kinetic energy for the C2 molecules to lie in the range of 2-9 eV. 

In a first set of experiments, each molecule is initially aligned orthogonal to the surface and shot along the direction 
of the [110] chains at an angle of 80° to the surface normal vector. The motivation of this particular choice was that 
the molecule may get adsorbed at the lowest energy position at the end of a C2n cluster, as identified in the diffusion 
studies. All the atoms, except the terminating hydrogen atoms on the bottom of the surface-slab were allowed to 
follow the Newtonian equations of motion. The molecules with up to 7 eV kinetic energy resulted in a C2 fragment 
on top of a C2n cluster, similar to the structure in Fig. ^ (b). At both 7 eV and 9 eV kinetic energy the molecule was 
deflected from the C2n cluster but was subsequently adsorbed as a lone C2 on a neighboring clean surface site. While 
the deflection at higher impact energy seems counter-intuitive at first sight, it must be recalled that the incident angle 
is high, so that the nature of the process is rather one of a steady dissipation of kinetic energy from the approaching 
C2 into the substrate until the molecule is slowed down enough to be deposited. 

In a second set of runs, a C2 was aligned parallel to the surface with an initial kinetic energy of 0.1-2 eV and 
a starting position on top of a Cg cluster. This leads again to the C2 fragment singly bonded to the Cg surface 
cluster, a configuration from which diffusion to either end is possible, as established before. There is a small region 
above the edge of a C2n (n=2,3) cluster, from where a C2, if given an initial velocity towards to the surface, can 
bond to the meta-stable minimum which precedes the diamond position, shown in Fig. ^. However, in the molecular 
dynamics simulations the barrier towards growth completion is too high in our time-scale, and we obtained solely the 
meta-stable minimum. 
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VII. SUMMARY AND CONCLUSIONS 



We have simulated diamond growth steps by subsequent deposition of C2 molecules onto a diamond (110) surface. 
We find that the initial C2 adsorption onto a clean (110) diamond surface proceeds with small barriers (0.1-0.2 eV) into 
the diamond lattice site. The growth mechanism and energetics of this insertion are similar to those on hydrogenated 
surfaces, as suggested previouslya Subsequent C2 additions on and around the initial adsorbate preferably leads to 
C2n chains forming along the [110] direction on the surface. The adsorption energies, as listed in Tab. || are in the 
range of 7-10 eV per C2 molecule at adsorption sites which lead to chain growth, and slightly smaller, 5-7 eV, for sites 
leading to defected growth. Some backbonds at the C2n chains are broken, leading to a graphene-like morphology, 
if a 50% coverage is reached for the added monolayer. However, the surface remains stable at this point. If the C2 
deposition continues, it induces healing of the broken backbonds due to re-formation of sp^ bonds at the terminus of 
the graphene sheets. 

We also find some meta-stable C2 defects during the growth, which may be responsible for starting new nucleation 
sites, a tendency that would explain the rather small grain size in the experimental studies which motivated this work. 
Low energy growth may be possible, if the approaching molecule has a kinetic energy within a window of 3-5 eV 
in order to overcome barriers and avoid defect formation. The direct adsorption into a diamond lattice position is 
possible only at the end of a C2n cluster or onto a clean site of the surface. Upon coalescence of different C2n chains, 
the remaining vacancies can be filled by the same growth species, although with slightly higher barriers (0.3-0.6 eV) 
than in the initial stages. 

We have also carried out surface diffusion studies. Because the C2 adsorptions normally result in tightly bonded 
adsorbate structures, inter-island diffusion of C2 molecule is rather unlikely. However, an intra-island diffusion path 
exists, where an added C2 molecule diffuses on top of a G2n chain until it reaches its end and is incorporated there. This 
diffusion behavior supports the C2 addition model which was evident from the deposition energetics. An implication 
of this fact is that an enhancement of surface diffusion rates would result in an increase in growth rate, by virtue of 
diffusion of migrating C2 species to and eventual incorporation into growth sites. 
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TABLE I. Calculated geometry of the clean relaxed C(llO) surface. Following the notation of Ref. dij is the bond 
length between atoms of the ith and jth layers; do is the equilibrium bulk bond length, and 9i is the bond angle within layer 
i. Ayi and Azi are the relaxation of the atomic positions in layer i along the y ([001], transversal to chains) and z ([110], 
surface-normal) directions, respectively. The units are A and degrees, respectively. 



Ref. 15 



this work 



do 

dii 

di2 

d22 
d23 

d-i3 



1.529 
1.419 
1.467 
1.490 
1.576 
1.526 



(100%) 
(-7.2%) 
(-4.1%) 
(-2.6%) 
(+3.1%) 
(-0.2%) 



1.544 
1.441 
1.490 
1.505 
1.587 
1.542 



(100%) 
(-6.7%) 
(-3.5%) 
(-2.6%) 
(+2.7%) 
(-0.2%) 



Azi 
Ay2 

AZ2 



123.3 
113.8 
±0.10 
-0.17 
±0.03 
+0.03 



122.0 
113.8 
±0.10 
-0.16 
±0.00 
+0.02 



TABLE IL Energy barriers (Ebarr) and adsorption energies (Eads) of C2 (and C) at the initial and final stages of growth, 
for varying target sites. The "top" position means that the initial position of the C2 is above a C^n cluster; "end" refers to an 
initial C2 position above the edge along the [110] direction of the C2n cluster, "same" and "other" express whether the added 
C2 is above the same or the adjacent (110) trough as the existing C2n cluster on the surface. Negative indices indicate missing 
C atoms in an otherwise contiguous chain or on the (110) surface. 
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C2 + 


(110):(2 X 1):C_2 


(110) 


(2x1) 


Fig. : 


(d) 


0.3 


-10.2" 




C2 + 


(110):(2 X 1) 


(110) 


(2 X 1):C2 


Fig. " 




0.0 


-7.2 




C2 + 


(110):C-3 


(110) 


C_i 






0.5 


-6.9 




C2 + 


(110):C-2 


(110) 








0.4 


-8.1 




C2 + 


(110):C_i 


(110):C 






0.0''/2.6 


-6.8*^7-8.3 




C + 


(110):C-i 


(110) 








0.0 


-10.4 





"C2n chain growth and coalescence processes 
''meta-stable state 
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TABLE III. The diffusion barriers Ebarr and the change in the total energy along the diffusion path. The low energy barriers 
(rows 4 and 5) are associated with the diffusion of a vertically aligned C2 on top of a C2n cluster. 



Path 


Ebarr /eV 


(Efl„ai-Ei„it)/eV 


C2 along valley 


3.8 


0.0 


C2 to other valley 


3.3 


0.0 


C2 to C4 along — » Cb 


3.7 


-1.8 


C2 to end of Ce — > Cs 


1.0"/0.6'' 


-3.7 


C2 to side of Ce — > Cs 


1.0^/0.6° 


-3.1 



''near center of Ce 

''near end of C2n 

'^from meta-stable to final energy minimum 
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FIG. 1. Overview of relaxed structures from subsequent depositions of C2 molecules onto a clean diamond (110) surface, 
(a) C2, (b) C4, (c) Ce, (d) Cg, all along the y = [001] direction, and (e) Cg, along the x = [110] direction. The numbers given 
are distances in A. Atoms of the adsorbate and the surface layer are shown bigger. Dark atoms indicate the adsorbate cluster, 
medium gray atoms its first and ring-forming second neighbors within the surface layer and lighter gray are other atoms. Small 
spheres at the bottom indicate hydrogen saturation. For (c)-(e), only partial models are shown. 




-8.1 eV 



FIG. 2. Initial steps for deposition of a C2 molecule onto a clean diamond (110) surface. The energies are given relative to 
a clean surface and a distant C2. Energies in parentheses indicate barriers. Atom designation is the same as in Fig. hi. 
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-7.2 eW 



-8.3 eV 



FIG. 3. Continued deposition of a C2 molecule onto a diamond (110) surface on sites next to an existing C2 adsorbate. The 
targeted neighboring site is (a) along the x = [llO] direction, and (b) along the y = [001] direction. The total energies are given 
relative to initially separated components. Energies in brackets indicate barriers. Atom designation is the same as in Fig. ^. 
Additional small markers indicate the target location. 




FIG. 4. Continued deposition of a C2 molecule onto a diamond (110) surface on top of an existing C2 adsorbate, with the 
transition state shown from the top and side. The resulting adsorbate cluster is topologically similar to a Ceo fragment. The 
panels show (a) the initial, (b) the transition, and (c) the final state. The total energies are given relative to initially separated 
components. Atom designation is the same as in Fig. 0. 
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FIG. 5. Continued deposition of a C2 molecule onto a diamond (110) surface over a Ce adsorbate with high insertion energy. 
The panels show (a) the initial, (b) the transition, and (c) the final state. The total energies are given relative to initially 
separated components. Atom designation is the same as in Fig. |l|, with the added C2 molecule shown in white. 



(a) C.3 



6.0 eV 




6.2 eV 




FIG. 6. Final stages of [110] chain growth and coalescence on a diamond (110) surface, (a) and (b) for a 3-site vacancy, (c) 
and (d) for a 2-site vacancy. Panels (a) and (c) show the meta-stable adsorption phases, and panels (b) and (d) the relaxed 
minimum configuration. Atom designation is the same as in Fig. |l|, with the added C2 molecule shown in white. The total 
energies are given relative to initially separated substrate and added C2 . Energies in parentheses indicate barriers. The dashed 
line in (b) indicates a cut used for the alternative view along [llO] in this panel. The arrow in (d) indicates an empirical bending 
radius for the graphene sheet. 
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FIG. 7. Graphitisation on a 50% covered diamond (110) surface in (2 x 1) reconstruction (a), and induced rebonding after 
deposition of C2, (b) in top view, (c) in side view. Dark, grey and white atoms indicate atoms in the top three monolayers, 
which are also shown in bigger size than the remaining atoms. The bottom monolayer is the hydrogen termination. 
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